Transcriptome profiles of Trypanosoma brucei rhodesiense in Malawi reveal focus specific gene expression profiles associated with pathology.

BACKGROUND: Sleeping sickness caused by Trypanosoma brucei rhodesiense is a fatal disease and endemic in Southern and Eastern Africa. There is an urgent need to develop novel diagnostic and control tools to achieve elimination of rhodesiense sleeping sickness which might be achieved through a better understanding of trypanosome gene expression and genetics using endemic isolates. Here, we describe transcriptome profiles and population structure of endemic T. b. rhodesiense isolates in human blood in Malawi. METHODOLOGY: Blood samples of r-HAT cases from Nkhotakota and Rumphi foci were collected in PaxGene tubes for RNA extraction before initiation of r-HAT treatment. 100 million reads were obtained per sample, reads were initially mapped to the human genome reference GRCh38 using HiSat2 and then the unmapped reads were mapped against Trypanosoma brucei reference transcriptome (TriTrypDB54_TbruceiTREU927) using HiSat2. Differential gene expression analysis was done using the DeSeq2 package in R. SNP calling from reads that were mapped to the T. brucei genome was done using GATK in order to identify T.b. rhodesiense population structure. RESULTS: 24 samples were collected from r-HAT cases of which 8 were from Rumphi and 16 from Nkhotakota foci. The isolates from Nkhotakota were enriched with transcripts for cell cycle arrest and stumpy form markers, whereas isolates in Rumphi focus were enriched with transcripts for folate biosynthesis and antigenic variation pathways. These parasite focus-specific transcriptome profiles are consistent with the more virulent disease observed in Rumphi and a less symptomatic disease in Nkhotakota associated with the non-dividing stumpy form. Interestingly, the Malawi T.b. rhodesiense isolates expressed genes enriched for reduced cell proliferation compared to the Uganda T.b. rhodesiense isolates. PCA analysis using SNPs called from the RNAseq data showed that T. b. rhodesiense parasites from Nkhotakota are genetically distinct from those collected in Rumphi. CONCLUSION: Our results suggest that the differences in disease presentation in the two foci is mainly driven by genetic differences in the parasites in the two major endemic foci of Rumphi and Nkhotakota rather than differences in the environment or host response.

Impact of environmental factors on Biomphalaria pfeifferi vector capacity leading to human infection by Schistosoma mansoni in two regions of western Côte d'Ivoire.

ABSTRACT: BACKGROUND: Intestinal schistosomiasis remains a worrying health problem, particularly in western Côte d'Ivoire, despite control efforts. It is therefore necessary to understand all the factors involved in the development of the disease, including biotic and abiotic factors. The aim of this study was to examine the factors that could support the maintenance of the intermediate host and its vectorial capacity in western Côte d'Ivoire. METHODS: Data on river physicochemical, microbiological, and climatic parameters, the presence or absence of snails with Schistosoma mansoni, and human infections were collected between January 2020 and February 2021. Spearman rank correlation tests, Mann-Whitney, analysis of variance (ANOVA), and an appropriate model selection procedure were used to analyze the data. RESULTS: The overall prevalence of infected snails was 56.05%, with infection reaching 100% in some collection sites and localities. Of 26 sites examined, 25 contained thermophilic coliforms and 22 contained Escherichia coli. Biomphalaria pfeifferi was observed in environments with lower land surface temperature (LST) and higher relative air humidity (RAH), and B. pfeifferi infection predominated in more acidic environments. Thermal coliforms and E. coli preferred higher pH levels. Lower maximum LST (LST_Max) and higher RAH and minimum LST (LST_Min) were favorable to E. coli, and lower LST_Max favored coliforms. The presence of B. pfeifferi was positively influenced by water temperature (T °C), LST_Min, RAH, and precipitation (Pp) (P < 0.05) and negatively influenced by pH, total dissolved solids (TDS), electrical conductivity (EC), LST_Max, and mean land surface temperature (LST). The parameters pH, TDS, EC, LST_Min, LST, and Pp had a positive impact on snail infection, while LST_Max had a negative impact on infection. Only pH had a positive effect on coliform and E. coli abundance. Of the 701 people examined for human schistosomiasis, 73.13% were positive for the point-of-care circulating cathodic antigen (POC-CCA) test and 12.01% for the Kato-Katz (KK) test. A positive correlation was established between human infections and the abundance of Biomphalaria (r2 = 0.879, P = 0.04959). CONCLUSIONS: The results obtained reflect the environmental conditions that are conducive to the maintenance of S. mansoni infection in this part of the country. To combat this infection as effectively as possible, it will be necessary not only to redouble efforts but also to prioritize control according to the level of endemicity at the village level.

Differences in gene expression profiles in early and late stage rhodesiense HAT individuals in Malawi.

T. b. rhodesiense is the causative agent of Rhodesian human African trypanosomiasis (r-HAT) in Malawi. Clinical presentation of r-HAT in Malawi varies between foci and differs from East African HAT clinical phenotypes. The purpose of this study was to gain more insights into the transcriptomic profiles of patients with early stage 1 and late stage 2 HAT disease in Malawi. Whole blood from individuals infected with T. b. rhodesiense was used for RNA-Seq. Control samples were from healthy trypanosome negative individuals matched on sex, age range, and disease foci. Illumina sequence FASTQ reads were aligned to the GRCh38 release 84 human genome sequence using HiSat2 and differential analysis was done in R Studio using the DESeq2 package. XGR, ExpressAnalyst and InnateDB algorithms were used for functional annotation and gene enrichment analysis of significant differentially expressed genes. RNA-seq was done on 23 r-HAT case samples and 28 healthy controls with 7 controls excluded for downstream analysis as outliers. A total of 4519 genes were significant differentially expressed (p adjusted <0.05) in individuals with early stage 1 r-HAT disease (n = 12) and 1824 genes in individuals with late stage 2 r-HAT disease (n = 11) compared to controls. Enrichment of innate immune response genes through neutrophil activation was identified in individuals with both early and late stages of the disease. Additionally, lipid metabolism genes were enriched in late stage 2 disease. We further identified uniquely upregulated genes (log2 Fold Change 1.4-2.0) in stage 1 (ZNF354C) and stage 2 (TCN1 and MAGI3) blood. Our data add to the current understanding of the human transcriptome profiles during T. b. rhodesiense infection. We further identified biological pathways and transcripts enriched than were enriched during stage 1 and stage 2 r-HAT. Lastly, we have identified transcripts which should be explored in future research whether they have potential of being used in combination with other markers for staging or r-HAT.

Transcriptome analysis of peripheral blood of Schistosoma mansoni infected children from the Albert Nile region in Uganda reveals genes implicated in fibrosis pathology.

Over 290 million people are infected by schistosomes worldwide. Schistosomiasis control efforts focus on mass drug treatment with praziquantel (PZQ), a drug that kills the adult worm of all Schistosoma species. Nonetheless, re-infections have continued to be detected in endemic areas with individuals living in the same area presenting with varying infection intensities. Our objective was to characterize the transcriptome profiles in peripheral blood of children between 10-15 years with varying intensities of Schistosoma mansoni infection living along the Albert Nile in Uganda. RNA extracted from peripheral blood collected from 44 S. mansoni infected (34 high and 10 low by circulating anodic antigen [CAA] level) and 20 uninfected children was sequenced using Illumina NovaSeq S4 and the reads aligned to the GRCh38 human genome. Differential gene expression analysis was done using DESeq2. Principal component analysis revealed clustering of gene expression by gender when S. mansoni infected children were compared with uninfected children. In addition, we identified 14 DEGs between S. mansoni infected and uninfected individuals, 56 DEGs between children with high infection intensity and uninfected individuals, 33 DEGs between those with high infection intensity and low infection intensity and no DEGs between those with low infection and uninfected individuals. We also observed upregulation and downregulation of some DEGs that are associated with fibrosis and its regulation. These data suggest expression of fibrosis associated genes as well as genes that regulate fibrosis in S. mansoni infection. The relatively few significant DEGS observed in children with schistosomiasis suggests that chronic S. mansoni infection is a stealth infection that does not stimulate a strong immune response.

Variants of IL6, IL10, FCN2, RNASE3, IL12B and IL17B loci are associated with Schistosoma mansoni worm burden in the Albert Nile region of Uganda.

BACKGROUND: Individuals genetically susceptible to high schistosomiasis worm burden may contribute disproportionately to transmission and could be prioritized for control. Identifying genes involved may guide development of therapy. METHODOLOGY/PRINCIPAL FINDINGS: A cohort of 606 children aged 10-15 years were recruited in the Albert Nile region of Uganda and assessed for Schistosoma mansoni worm burden using the Up-Converting Particle Lateral Flow (UCP-LF) test detecting circulating anodic antigen (CAA), point-of-care Circulating Cathodic Antigen (POC-CCA) and Kato-Katz tests. Whole genome genotyping was conducted on 326 children comprising the top and bottom 25% of worm burden. Linear models were fitted to identify variants associated with worm burden in preselected candidate genes. Expression quantitative trait locus (eQTL) analysis was conducted for candidate genes with UCP-LF worm burden included as a covariate. Single Nucleotide Polymorphism loci associated with UCP-LF CAA included IL6 rs2066992 (OR = 0.43, p = 0.0006) and rs7793163 (OR = 2.0, p = 0.0007); IL21 SNP kgp513476 (OR 1.79, p = 0.0025) and IL17B SNP kgp708159 (OR = 0.35, p = 0.0028). A haplotype in the IL10 locus was associated with lower worm burden (OR = 0.53, p = 0.015) and overlapped SNPs rs1800896, rs1800871 and rs1800872. Significant haplotypes (p<0.05, overlapping significant SNP) associated with worm burden were observed in IL6 and the Th17 pathway IL12B and IL17B genes. There were significant eQTL in the IL6, IL5, IL21, IL25 and IFNG regions. CONCLUSIONS: Variants associated with S. mansoni worm burden were in IL6, FCN2, RNASE3, IL10, IL12B and IL17B gene loci. However only eQTL associations remained significant after Bonferroni correction. In summary, immune balance, pathogen recognition and Th17 pathways may play a role in modulating Schistosoma worm burden. Individuals carrying risk variants may be targeted first in allocation of control efforts to reduce the burden of schistosomiasis in the community.

Association between polymorphisms of IL4, IL13, IL10, STAT6 and IFNG genes, cytokines and immunoglobulin E levels with high burden of Schistosoma mansoni in children from schistosomiasis endemic areas of Cameroon.

Eliminating schistosomiasis as a public health problem by 2030 requires a better understanding of the disease transmission, especially the asymmetric distribution of worm burden in individuals living and sharing the same environment. It is in this light that this study was designed to identify human genetic determinants associated with high burden of S. mansoni and also with the plasma concentrations of IgE and four cytokines in children from two schistosomiasis endemic areas of Cameroon. In school-aged children of schistosomiasis endemic areas of Makenene and Nom-Kandi of Cameroon, S. mansoni infections and their infection intensities were evaluated in urine and stool samples using respectively the Point-of-care Circulating Cathodic Antigen test (POC-CCA) and the Kato Katz (KK) test. Thereafter, blood samples were collected in children harbouring high burden of schistosome infections as well as in their parents and siblings. DNA extracts and plasma were obtained from blood. Polymorphisms at 14 loci of five genes were assessed using PCR-restriction fragment length polymorphism and amplification-refractory mutation system. The ELISA test enabled to determine the plasma concentrations of IgE, IL-13, IL-10, IL-4 and IFN-γ. The prevalence of S. mansoni infections was significantly higher (P < 0.0001 for POC-CCA; P = 0.001 for KK) in Makenene (48.6% for POC-CCA and 7.9% for KK) compared to Nom-Kandi (31% for POC-CCA and 4.3% for KK). The infection intensities were also higher (P < 0.0001 for POC-CCA; P = 0.001 for KK) in children from Makenene than those from Nom-Kandi. The allele C of SNP rs3024974 of STAT6 was associated with an increased risk of bearing high burden of S. mansoni both in the additive (p = 0.009) and recessive model (p = 0.01) while the allele C of SNP rs1800871 of IL10 was protective (p = 0.0009) against high burden of S. mansoni. The alleles A of SNP rs2069739 of IL13 and G of SNP rs2243283 of IL4 were associated with an increased risk of having low plasma concentrations of IL-13 (P = 0.04) and IL-10 (P = 0.04), respectively. This study showed that host genetic polymorphisms may influence the outcome (high or low worm burden) of S. mansoni infections and also the plasma concentrations of some cytokines.

High prevalence of Schistosoma mansoni infection and stunting among school age children in communities along the Albert-Nile, Northern Uganda: A cross sectional study.

BACKGROUND: Knowing the prevalence of schistosomiasis is key to informing programmes to control and eliminate the disease as a public health problem. It is also important to understand the impact of infection on child growth and development in order to allocate appropriate resources and effort to the control of the disease. METHODS: We conducted a survey to estimate the prevalence of schistosomiasis among school aged children in villages along the Albert-Nile shore line in the district of Pakwach, North Western Uganda. A total of 914 children aged between 10-15 years were screened for Schistosoma mansoni using the POC-CCA and Kato Katz (KK) techniques. The infection intensities were assessed by POC-CCA and KK as well as CAA tests. The KK intensities were also correlated with POC-CCA and with CAA intensity. Anthropometric measurements were also taken and multivariate analysis was carried out to investigate their association with infection status. RESULTS: The prevalence of schistosomiasis using the POC-CCA diagnostic test was estimated at 85% (95% CI: 83-87), being highest amongst children living closer to the Albert-Nile shoreline. Visual scoring of the POC-CCA results was more sensitive than the Kato Katz test and was positively correlated with the quantified infection intensities by the CAA test. The majority of the children were underweight (BMI<18.5), and most notably, boys had significantly lower height for age (stunting) than girls in the same age range (p < 0.0001), but this was not directly associated with S. mansoni infection. CONCLUSION: High prevalence of S. mansoni infection in the region calls for more frequent mass drug administration with praziquantel. We observed high levels of stunting which was not associated with schistosomiasis. There is a need for improved nutrition among the children in the area.

Field assessment in Cameroon of a reader of POC-CCA lateral flow strips for the quantification of Schistosoma mansoni circulating cathodic antigen in urine.

BACKGROUND: Determining Schistosoma mansoni infection rate and intensity is challenging due to the low sensitivity of the Kato-Katz (KK) test that underestimates the true disease prevalence. Circulating cathodic antigen (CCA) excreted in urine is constantly produced by adult worms and has been used as the basis of a simple, non-invasive point of care test (POC-CCA) for Schistosoma mansoni infections. Although the abundance of CCA in urine is proportional to worm burden, the POC-CCA test is marketed as a qualitative test, making it difficult to investigate the wide range of infection intensities. This study was designed to compare the prevalence and intensity of S. mansoni by KK and POC-CCA and quantify, on fresh and frozen (<-20°C) urine samples, CCA using the visual scores and the ESEquant LR3 reader. METHODOLOGY: Stool and urine samples were collected from 759 school-aged children. The prevalence and intensity of S. mansoni were determined using KK and POC-CCA. The degree of the positivity of POC-CCA was estimated by quantifying CCA on fresh and frozen urine samples using visual scores and strip reader. The prevalence, the infection intensity as well the relative amounts of CCA were compared. RESULTS: The S. mansoni infection rates inferred from POC-CCA and KK were 40.7% and 9.4% respectively. Good correlations were observed between infection intensities recorded by; i) the reader and visual scoring system on fresh (Rho = 0.89) and frozen samples (Rho = 0.97), ii) the reader on fresh urine samples and KK (epg) (Rho = 0.44). Nevertheless, 238 POC-CCA positive children were negative for KK, and sixteen of them had high levels of CCA. The correlation between results from the reader on fresh and frozen samples was good (Rho = 0.85). On frozen samples, CCA was not detected in 55 samples that were positive in fresh urine samples. CONCLUSION: This study confirmed the low sensitivity of KK and the high capacity of POC-CCA to provide reliable data on the prevalence and intensity of S. mansoni infections. The lateral flow reader enabled accurate quantification of CCA under field conditions on fresh and frozen urine samples with less time and effort than KK.

The Genetics of Human Schistosomiasis Infection Intensity and Liver Disease: A Review.

Schistosomiasis remains the fourth most prevalent parasitic disease affecting over 200 million people worldwide. Control efforts have focussed on the disruption of the life cycle targeting the parasite, vector and human host. Parasite burdens are highly skewed, and the majority of eggs are shed into the environment by a minority of the infected population. Most morbidity results from hepatic fibrosis leading to portal hypertension and is not well-correlated with worm burden. Genetics as well as environmental factors may play a role in these skewed distributions and understanding the genetic risk factors for intensity of infection and morbidity may help improve control measures. In this review, we focus on how genetic factors may influence parasite load, hepatic fibrosis and portal hypertension. We found 28 studies on the genetics of human infection and 20 studies on the genetics of pathology in humans. S. mansoni and S. haematobium infection intensity have been showed to be controlled by a major quantitative trait locus SM1, on chromosome 5q31-q33 containing several genes involved in the Th2 immune response, and three other loci of smaller effect on chromosomes 1, 6, and 7. The most common pathology associated with schistosomiasis is hepatic and portal vein fibroses and the SM2 quantitative trait locus on chromosome six has been linked to intensity of fibrosis. Although there has been an emphasis on Th2 cytokines in candidate gene studies, we found that four of the five QTL regions contain Th17 pathway genes that have been included in schistosomiasis studies: IL17B and IL12B in SM1, IL17A and IL17F in 6p21-q2, IL6R in 1p21-q23 and IL22RA2 in SM2. The Th17 pathway is known to be involved in response to schistosome infection and hepatic fibrosis but variants in this pathway have not been tested for any effect on the regulation of these phenotypes. These should be priorities for future studies.

Candidate gene family-based and case-control studies of susceptibility to high Schistosoma mansoni worm burden in African children: a protocol.

Background: Approximately 25% of the risk of Schistosoma mansoni is associated with host genetic variation. We will test 24 candidate genes, mainly in the T h2 and T h17 pathways, for association with S. mansoni infection intensity in four African countries, using family based and case-control approaches. Methods: Children aged 5-15 years will be recruited in S. mansoni endemic areas of Ivory Coast, Cameroon, Uganda and the Democratic Republic of Congo (DRC). We will use family based (study 1) and case-control (study 2) designs. Study 1 will take place in Ivory Coast, Cameroon, Uganda and the DRC. We aim to recruit 100 high worm burden families from each country except Uganda, where a previous study recruited at least 40 families. For phenotyping, cases will be defined as the 20% of children in each community with heaviest worm burdens as measured by the circulating cathodic antigen (CCA) assay. Study 2 will take place in Uganda. We will recruit 500 children in a highly endemic community. For phenotyping, cases will be defined as the 20% of children with heaviest worm burdens as measured by the CAA assay, while controls will be the 20% of infected children with the lightest worm burdens. Deoxyribonucleic acid (DNA) will be genotyped on the Illumina H3Africa SNP (single nucleotide polymorphisms) chip and genotypes will be converted to sets of haplotypes that span the gene region for analysis. We have selected 24 genes for genotyping that are mainly in the Th2 and Th17 pathways and that have variants that have been demonstrated to be or could be associated with Schistosoma infection intensity.   Analysis: In the family-based design, we will identify SNP haplotypes disproportionately transmitted to children with high worm burden. Case-control analysis will detect overrepresentation of haplotypes in extreme phenotypes with correction for relatedness by using whole genome principal components.

High-depth African genomes inform human migration and health.

The African continent is regarded as the cradle of modern humans and African genomes contain more genetic variation than those from any other continent, yet only a fraction of the genetic diversity among African individuals has been surveyed1. Here we performed whole-genome sequencing analyses of 426 individuals-comprising 50 ethnolinguistic groups, including previously unsampled populations-to explore the breadth of genomic diversity across Africa. We uncovered more than 3 million previously undescribed variants, most of which were found among individuals from newly sampled ethnolinguistic groups, as well as 62 previously unreported loci that are under strong selection, which were predominantly found in genes that are involved in viral immunity, DNA repair and metabolism. We observed complex patterns of ancestral admixture and putative-damaging and novel variation, both within and between populations, alongside evidence that Zambia was a likely intermediate site along the routes of expansion of Bantu-speaking populations. Pathogenic variants in genes that are currently characterized as medically relevant were uncommon-but in other genes, variants denoted as 'likely pathogenic' in the ClinVar database were commonly observed. Collectively, these findings refine our current understanding of continental migration, identify gene flow and the response to human disease as strong drivers of genome-level population variation, and underscore the scientific imperative for a broader characterization of the genomic diversity of African individuals to understand human ancestry and improve health.

High Levels of Genetic Diversity within Nilo-Saharan Populations: Implications for Human Adaptation.

Africa contains more human genetic variation than any other continent, but the majority of the population-scale analyses of the African peoples have focused on just two of the four major linguistic groups, the Niger-Congo and Afro-Asiatic, leaving the Nilo-Saharan and Khoisan populations under-represented. In order to assess genetic variation and signatures of selection within a Nilo-Saharan population and between the Nilo-Saharan and Niger-Congo and Afro-Asiatic, we sequenced 50 genomes from the Nilo-Saharan Lugbara population of North-West Uganda and 250 genomes from 6 previously unsequenced Niger-Congo populations. We compared these data to data from a further 16 Eurasian and African populations including the Gumuz, another putative Nilo-Saharan population from Ethiopia. Of the 21 million variants identified in the Nilo-Saharan population, 3.57 million (17%) were not represented in dbSNP and included predicted non-synonymous mutations with possible phenotypic effects. We found greater genetic differentiation between the Nilo-Saharan Lugbara and Gumuz populations than between any two Afro-Asiatic or Niger-Congo populations. F3 tests showed that Gumuz contributed a genetic component to most Niger-Congo B populations whereas Lugabara did not. We scanned the genomes of the Lugbara for evidence of selective sweeps. We found selective sweeps at four loci (SLC24A5, SNX13, TYRP1, and UVRAG) associated with skin pigmentation, three of which already have been reported to be under selection. These selective sweeps point toward adaptations to the intense UV radiation of the Sahel.

Copy number variation in human genomes from three major ethno-linguistic groups in Africa.

BACKGROUND: Copy number variation is an important class of genomic variation that has been reported in 75% of the human genome. However, it is underreported in African populations. Copy number variants (CNVs) could have important impacts on disease susceptibility and environmental adaptation. To describe CNVs and their possible impacts in Africans, we sequenced genomes of 232 individuals from three major African ethno-linguistic groups: (1) Niger Congo A from Guinea and Côte d'Ivoire, (2) Niger Congo B from Uganda and the Democratic Republic of Congo and (3) Nilo-Saharans from Uganda. We used GenomeSTRiP and cn.MOPS to identify copy number variant regions (CNVRs). RESULTS: We detected 7608 CNVRs, of which 2172 were only deletions, 2384 were only insertions and 3052 had both. We detected 224 previously un-described CNVRs. The majority of novel CNVRs were present at low frequency and were not shared between populations. We tested for evidence of selection associated with CNVs and also for population structure. Signatures of selection identified previously, using SNPs from the same populations, were overrepresented in CNVRs. When CNVs were tagged with SNP haplotypes to identify SNPs that could predict the presence of CNVs, we identified haplotypes tagging 3096 CNVRs, 372 CNVRs had SNPs with evidence of selection (iHS > 3) and 222 CNVRs had both. This was more than expected (p < 0.0001) and included loci where CNVs have previously been associated with HIV, Rhesus D and preeclampsia. When integrated with 1000 Genomes CNV data, we replicated their observation of population stratification by continent but no clustering by populations within Africa, despite inclusion of Nilo-Saharans and Niger-Congo populations within our dataset. CONCLUSIONS: Novel CNVRs in the current study increase representation of African diversity in the database of genomic variants. Over-representation of CNVRs in SNP signatures of selection and an excess of SNPs that both tag CNVs and are subject to selection show that CNVs may be the actual targets of selection at some loci. However, unlike SNPs, CNVs alone do not resolve African ethno-linguistic groups. Tag haplotypes for CNVs identified may be useful in predicting African CNVs in future studies where only SNP data is available.

The distribution of haemoglobin C and its prevalence in newborns in Africa.

Haemoglobin C (HbC) is one of the commonest structural haemoglobin variants in human populations. Although HbC causes mild clinical complications, its diagnosis and genetic counselling are important to prevent inheritance with other haemoglobinopathies. Little is known about its contemporary distribution and the number of newborns affected. We assembled a global database of population surveys. We then used a Bayesian geostatistical model to create maps of HbC frequency across Africa and paired our predictions with high-resolution demographics to calculate heterozygous (AC) and homozygous (CC) newborn estimates and their associated uncertainty. Data were too sparse outside Africa for this methodology to be applied. The highest frequencies were found in West Africa but HbC was commonly found in other parts of the continent. The expected annual numbers of AC and CC newborns in Africa were 672,117 (interquartile range (IQR): 642,116-705,163) and 28,703 (IQR: 26,027-31,958), respectively. These numbers are about two times previous estimates.

Online biomedical resources for malaria-related red cell disorders.

Warnings about the expected increase of the global public health burden of malaria-related red cell disorders are accruing. Past and present epidemiological data are necessary to track spatial and temporal changes in the frequencies of these genetic disorders. A number of open access biomedical databases including data on malaria-related red cell disorders have been launched over the last two decades. Here, we review the content of these databases, most of which focus on genetic diversity, and we describe a new epidemiological resource developed by the Malaria Atlas Project. To tackle upcoming public health challenges, the integration of epidemiological and genetic data is important. As many countries are considering implementing national screening programs, strategies to make such data more accessible are also needed.

Global epidemiology of sickle haemoglobin in neonates: a contemporary geostatistical model-based map and population estimates.

BACKGROUND: Reliable estimates of populations affected by diseases are necessary to guide efficient allocation of public health resources. Sickle haemoglobin (HbS) is the most common and clinically significant haemoglobin structural variant, but no contemporary estimates exist of the global populations affected. Moreover, the precision of available national estimates of heterozygous (AS) and homozygous (SS) neonates is unknown. We aimed to provide evidence-based estimates at various scales, with uncertainty measures. METHODS: Using a database of sickle haemoglobin surveys, we created a contemporary global map of HbS allele frequency distribution within a Bayesian geostatistical model. The pairing of this map with demographic data enabled calculation of global, regional, and national estimates of the annual number of AS and SS neonates. Subnational estimates were also calculated in data-rich areas. FINDINGS: Our map shows subnational spatial heterogeneities and high allele frequencies across most of sub-Saharan Africa, the Middle East, and India, as well as gene flow following migrations to western Europe and the eastern coast of the Americas. Accounting for local heterogeneities and demographic factors, we estimated that the global number of neonates affected by HbS in 2010 included 5,476,000 (IQR 5,291,000-5,679,000) AS neonates and 312,000 (294,000-330,000) SS neonates. These global estimates are higher than previous conservative estimates. Important differences predicted at the national level are discussed. INTERPRETATION: HbS will have an increasing effect on public health systems. Our estimates can help countries and the international community gauge the need for appropriate diagnoses and genetic counselling to reduce the number of neonates affected. Similar mapping and modelling methods could be used for other inherited disorders. FUNDING: The Wellcome Trust.

G6PD deficiency prevalence and estimates of affected populations in malaria endemic countries: a geostatistical model-based map.

BACKGROUND: Primaquine is a key drug for malaria elimination. In addition to being the only drug active against the dormant relapsing forms of Plasmodium vivax, primaquine is the sole effective treatment of infectious P. falciparum gametocytes, and may interrupt transmission and help contain the spread of artemisinin resistance. However, primaquine can trigger haemolysis in patients with a deficiency in glucose-6-phosphate dehydrogenase (G6PDd). Poor information is available about the distribution of individuals at risk of primaquine-induced haemolysis. We present a continuous evidence-based prevalence map of G6PDd and estimates of affected populations, together with a national index of relative haemolytic risk. METHODS AND FINDINGS: Representative community surveys of phenotypic G6PDd prevalence were identified for 1,734 spatially unique sites. These surveys formed the evidence-base for a Bayesian geostatistical model adapted to the gene's X-linked inheritance, which predicted a G6PDd allele frequency map across malaria endemic countries (MECs) and generated population-weighted estimates of affected populations. Highest median prevalence (peaking at 32.5%) was predicted across sub-Saharan Africa and the Arabian Peninsula. Although G6PDd prevalence was generally lower across central and southeast Asia, rarely exceeding 20%, the majority of G6PDd individuals (67.5% median estimate) were from Asian countries. We estimated a G6PDd allele frequency of 8.0% (interquartile range: 7.4-8.8) across MECs, and 5.3% (4.4-6.7) within malaria-eliminating countries. The reliability of the map is contingent on the underlying data informing the model; population heterogeneity can only be represented by the available surveys, and important weaknesses exist in the map across data-sparse regions. Uncertainty metrics are used to quantify some aspects of these limitations in the map. Finally, we assembled a database of G6PDd variant occurrences to inform a national-level index of relative G6PDd haemolytic risk. Asian countries, where variants were most severe, had the highest relative risks from G6PDd. CONCLUSIONS: G6PDd is widespread and spatially heterogeneous across most MECs where primaquine would be valuable for malaria control and elimination. The maps and population estimates presented here reflect potential risk of primaquine-associated harm. In the absence of non-toxic alternatives to primaquine, these results represent additional evidence to help inform safe use of this valuable, yet dangerous, component of the malaria-elimination toolkit. Please see later in the article for the Editors' Summary.

The global distribution of the Duffy blood group.

Blood group variants are characteristic of population groups, and can show conspicuous geographic patterns. Interest in the global prevalence of the Duffy blood group variants is multidisciplinary, but of particular importance to malariologists due to the resistance generally conferred by the Duffy-negative phenotype against Plasmodium vivax infection. Here we collate an extensive geo-database of surveys, forming the evidence-base for a multi-locus Bayesian geostatistical model to generate global frequency maps of the common Duffy alleles to refine the global cartography of the common Duffy variants. We show that the most prevalent allele globally was FY*A, while across sub-Saharan Africa the predominant allele was the silent FY*B(ES) variant, commonly reaching fixation across stretches of the continent. The maps presented not only represent the first spatially and genetically comprehensive description of variation at this locus, but also constitute an advance towards understanding the transmission patterns of the neglected P. vivax malaria parasite.

Global distribution of the sickle cell gene and geographical confirmation of the malaria hypothesis.

It has been 100 years since the first report of sickle haemoglobin (HbS). More than 50 years ago, it was suggested that the gene responsible for this disorder could reach high frequencies because of resistance conferred against malaria by the heterozygous carrier state. This traditional example of balancing selection is known as the 'malaria hypothesis'. However, the geographical relationship between the transmission intensity of malaria and associated HbS burden has never been formally investigated on a global scale. Here, we use a comprehensive data assembly of HbS allele frequencies to generate the first evidence-based map of the worldwide distribution of the gene in a Bayesian geostatistical framework. We compare this map with the pre-intervention distribution of malaria endemicity, using a novel geostatistical area-mean comparison. We find geographical support for the malaria hypothesis globally; the relationship is relatively strong in Africa but cannot be resolved in the Americas or in Asia.